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ABSTRACT 

We present a detailed investigation of the near-to-far infrared (IR) spectral energy distribution 
(SED) of a large sample of Spitzer-observed active galactic nuclei (AGN). We fitted the spec- 
tra of 51 narrow line Seyfert 1 galaxies (NLSls) and 64 broad line Seyfert 1 galaxies (BLSls) 
using a three component model: a clumpy torus, a dusty narrow line region (NLR) and hot 
pure-graphite dust located in the outer part of the broad line region (BLR). The fitting is per- 
formed on star formation (SF) subtracted SEDs using SF templates that take into account the 
entire range of possible host galaxy properties. We find that the mid-IR intrinsic emission of 
NLSls and BLSls are very similar, regardless of the AGN luminosity, with long wavelength 
downturn at around 20-25 jUm. We present a detailed model of the hot dust component that 
takes into account the distribution of dust temperature within the clouds and their emission 
line spectrum. The hot dust continuum provides a very good fit to the observed near-IR con- 
tinuum spectrum. Most line emission in this component is dramatically suppressed, except 
Mg II A2798 and He 1 lines that are still contributing significantly to the total BLR spectrum. 
We calculate the covering factors of all the AGN components and show that the covering fac- 
tor of the hot-dust clouds is about 0.15-0.35, similar to the covering factor of the torus, and is 
anti-correlated with the source luminosity and the normalized accretion rate. 
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1 INTRODUCTION 

The unification scheme of active galactic nuclei (AGN) requires an 
anisotropic obscuring structure that surrounds the central accret- 



ing black hole (e.g.. iKrolik & BegelmaiJl 19881 ; lAntonuccilll993h . 
In this picture, the bulk of the radiation from the central engine 
is absorbed by the obscuring structure and re-emitted mainly at 
mid-infrared (MIR) wavelengths. The infrared (IR) spectral en- 
ergy distribution (SED) of AGN also shows emission by star form- 
ing (SF) regions in the host galaxy. Recent studies have shown 
that polycyclic aromatic hydrocarbon (PA H) features can serve as 
good indicators for this SF activity (e.g., ISchweitzer et aljflOOd : 
iLaMassa et"^l2010l) . Measuring and subtracting the SF contribu- 
tion to the IR SED enables to focus on the AGN-produced SED. 
iNetzeretalJbOOTh used the PAH features to subtract the SF contri- 
bution from the average SED of the QUEST sample of PG quasars 
(QSOs). This was done for two sub groups of low and high far-IR 
(FIR) luminosity. The two subtracted SEDs were found to be re- 
markabl y similar suggesting that most AGN have similar intrinsic 
SEDs. lDeo et al.l 1200 9) reached a similar conclusion while inves- 
tigating samples of both type-I and type-II Seyfert galaxies. Both 
studies found that the intrinsic AGN SEDs exhibit two main bumps 
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around ^ 10 and 20-/im. These are associated with broad sili- 
cate emission features as well as emi ssion from cooler dust clouds. 
More recently Mulla nev et al.l ( 1201 ih found that the second bump 
may peak at somewhat longer wavelengths (~ 30-/xm) indicating 
the presence of even colder dust. 

A main component of the obscuring circumnuclear struc- 
ture is believed to be a dusty torus. The MIR SED of such a 
torus depends on its dimensions and geometry, the density dis- 
tribution and the dust grain properties. Several attempts have 
been made t o model such tori assuming smooth density distri 
butions (e.g ., 'Pier & Krolik"l992','l993';'Granato & Danese 199' 
Efstathiou & Rowan-Robinson 1995; van Bemmel & Dullemonj 
20031 ; Schartmann et al. 2005) and clumpy density distributions 
(e.g..'"Nen kova et al..2008a ; .Honig et al.l2010l) . Other AGN compo- 
nents can contribute to the observed MIR spectrum. Some of this 
emission may originate farther from the central radiation source, 
at distances exceeding the dimensions of the torus. Dusty gas in 
the narrow-line region (NLR) may be the source of such radiation 
jSchweitzer etai]|2008l ; iMor et alj[2009l. hereafter M09). Thus, a 
significant contribution at ^ 2Q-jxm due to components not related 
to the torus, must be considered. 

Dust reverberation measurements of several nearby AGN, 
based on the V and K band emission, lead to the conclusion that 
the near IR (NIR) emission in these sources is dominated by ther- 
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mal radiation from hot dust v ery close to the centre jMinezaki et alJ 
l2004ISuganuma et al.ll2006h . There have been several studies that 
fitted the NIR to MIR SEDs of AGN using a black- body spectrum 
to represent emission f r om such hot dust (e.g., Ede lson & MalkanI 

, Ri ffeletal.ll200a 
etalj j20Ilh found 



Il98a lBarvaim^ ll987l ; iKishimoto et"aL i]2007: 
M09; iDeo et alj|201ll) . More recently, tandt ( 
similar results by fitting the optical to NIR SED of 23 AGN. The 
modelled temperature of this component, in all these studies, is 
high ( ^ 1200 K), regardless of the AGN lumi nosity, and is consis- 
tent w ith pure-graphite dust emission (M09). iMor & TrakhtenbrotI 
( l201lh used the (observed) optical to MIR SEDs of ~ 15000 type-I 
AGN to conclude that a hot pure-graphite dust component is ob- 
served in the vast majority of such sources. Several studies found 
that the NIR lu minosity strongly correlates with the luminosity o f 
the AGN (e.g., iGallagher et al.ll2007l ; iMor & TrakhtenbrotI 1201 ll) . 
[Gallagher et al.l ( 2007 ) concluded that this tight relation suggests a 
consta nt covering factor of dusty clouds regardless of AGN lumi- 
nosity. iMor & TrakhtenbrotI ( 1201 ih found that the covering factor 
of the hot dust component decreases with increasing AGN bolo- 
metric luminosity. 

Following M09, we aim to investigate a larger sample of type-I 
AGN with a broad range of luminosities, black hole masses (Mge), 
and normalized accretion rates {L/L^^i). We fit the observed MIR 
spectra of 115 type-I AGN using a 3-component model made of 
a clumpy torus, dusty NLR clouds and very hot dust clouds. Our 
main goal is to explore the physical properties of these components, 
focusing mainly on the hot dust. The high quality spectra made 
available by the IRS spectrom eter on-board the Spitzer Space Ob- 
servatory l lHouck et alj|2004l) allows us to explore these ideas. In 
^we describe the observational data. In ^we detail our model 
and the fitting procedure and in ^we present the results of this 
procedure and discuss their implications. 



2 SAMPLE SELECTION OBSERVATIONS AND DATA 
REDUCTION 

Our database consists of local type-I AGN collected from several 
samples, covering a wide range in luminosity and divided into two 
groups of NLS Is and BLS Is, according to their FWHM[H/3], with 
a division line at 2000 km s"'. 20 NLSls in our sample are taken 
from the ROSAT sample of nearby AGN ( IThomas etal]|l998h . and 
are a part of the Spit zer PID 20241 (PI D. Lutz) recently described 
by ISani et al. hereafter SIO). SIO supplemented these ob- 

jects by a large number of additional SpitzerflRS archival obser- 
vations of NLSls and BLSls. These additional objects were se- 
lected from the 12th edit ion of the Catalog of Quasar s and Ac- 
tive Nuclei compiled by IVeron- Cettv & Veroril ( l2006l) . We also 
added to the sample all AGN in the QUEST Spitzer spectroscopy 
project (PID 3187, PI S . Veilleux), which is described in detail in 
ISchweitzer et alj ( l2006h . iNetzer et all ( |2007|) an dM09. Most of the 
QUE ST objects are Palomar- Green (PG) QSO s dSchmidt & GreerJ 
Il983l) and are taken from Guv on etal] j2006l) . In total, our sample 
consists of 1 15 sources (51 NLSls, 64 BLSls) spanning a luminos- 



ity range of L5io()~ lO'*^'^^'*^'^ ergs^' where L5100 stands for XL^ 
at rest wavelength 5100 A. 

The Spitzer observations and the data reduction proce- 
dure of the SIO and Q UEST samples are detailed in SIO and 
ISchweitzer et alj ( |2006[) . respectively and are summarized here for 
completion. IRS spectra for all objects were taken in either low or 
high resolution modes and cover the wavelength range of 5-35 flm 
in the observed frame. The standard slit widths of 3". 6 to 1 l".l in- 



clude flux from the host galaxies and the vicinity of the AGN. We 
started the data reduction from the basic calibrated data (BCD) pro- 
vided by the Spitzer pipeline. Specially developed IDL-based tools 
are then used for removing outlying val ues for individual pi xels and 
for sky subtraction. The SMART tool ( iHigdon et al]|2004t) is used 
for extraction of the final spectra. 



We supplemented the Spitzer spectra with NIR data ob- 
tained from the 2MA SS extended and point source catalogues 
jSkrutskie et al.l2006h . Approximately 2/3 of the sources have been 
detected by the 2MASS survey. We re-measured the 2MASS im- 
ages using an IRS-slit-like aperture in order to avoid biases regard- 
ing the large aperture used in the 2MASS catalogues. This gives a 
better calibration to the NIR measurements with respect to the spec- 
trum of the source, and is more relevant for the extended sources in 
our sample ( 50% of the sample). The main source of uncertainty 
related to the NIR data is due to the fact that many of these images 
were taken more than 20 years before the Spitzer observations. For 
the longer wavelengths of th e SED we use IRAS 60-/im photome- 
try jNeugebauer et al.lll984h . Aperture corrections, similar to those 
done for the 2MASS data, are not possible for the IRAS data due to 
the low spatial resolution of the IRAS instrument (~ 2'). This con- 
tributes to the uncertainty in the FIR flux of the extended objects in 
our sample and will be discussed in 33.11 



In order to compute L[,oi for all sources we adopt 
the iMarconi et al. I I2OO4) ' 'intrinsic", luminosity dependant 
SED, which provides a polynomial prescription for estimating 
AL;^(4400-A) at every Lboj. We use this prescription for L5100 
by adopting a /y °= v^"^ power-law approximation. For the SIO 
sample, values of Lsioo were obtained from various works (see 
SIO and references therein) For most objects in the QUEST sam- 
ple the values of L-jinn are obtained from the observations of 
iBoroson & Greed jl992t) . For PG 1244-H026, PG 1001-1-054, and 
PG 01 57-1-001 we used spectra from the s eventh data release of 
the Sloan Digital Sky Survey (SDSS/DR7,'Abazaiian et al."2009h 
that w ere measured in the way described in Netzer & Trakhtenbroj 
j2007h . We have estimated Mbh and L/LEdd for all sources using 
the procedure described in Netzer & Trakhtenbrot (2007). In this 
procedure (the "virial" mass determination, see Eq. 1 in Netzer & 
Trakhtenbrot 2007), L5100 and FWHM(Hj8) are combined to obtain 
Mbh- I^/^Edd is obtained by using the adopted bolometric correc- 
tion factor (BC). 



There are two uncertainties associated with the use of L(,oi. 
The first is source variability, which is an important effect since the 
optical and MIR observations were separated by many years. We 
estimate this uncertainty to be a factor of ~ 1.5. The second un- 
certainty involves the approximation used for the BC. We estimate 
this uncertainty to be ~30%. Both uncertainties affect the derived 
model parameters such as the covering factors and the distances to 
the NLR and hot dust clouds. Because of this, we do not attach great 
importance to specific values obtained from our best acceptable 
models. On the other hand, the sample is large enough to enable a 
significant analysis of its mean properties since the larger of these 
effects, due to source variability, is likely to be random. The uncer- 
tainty due to the estimate of BC is smaller but more problematic 
since the expression we use can under-estimate or over-estimate 
Lbol for all sources. This can introduce systematic differences, e.g. 
a smaller covering factors for all sources. 
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Figure 1. Ratio between FIR and MIR fluxes vs. IR luminosity. Gray points 
represent IR bright galaxies from the The IRAS Revised Bright Galaxy 
Sample ISan ders et al. (2003). Calculated ratios from theoretical SF tem- 
plates of lCharv & Elbazl (2001) are presented as triangles and the observed 
ratio in the 5 galaxies used as SF templates in this paper are shown as black 
squares. The horizontal dashed lines show flux ratios of 3 and 13. 



3 SPECTRAL DECOMPOSITION 

One of the main goals of the present work is to investigate the 
2-60 flm SED of the sources in our sample, following the gen- 
eral scheme of M09. Our models include different physical com- 
ponents: a very hot pure-graphite dust, a dusty torus containing 
graphite and silicate grains, a dusty NLR, and a SF host galaxy. 
We start by separating the AGN-related components from the host 
galaxy component. The subtraction of the SF contribution is a ma- 
jor challenge and a major source of uncertainty mainly at longer 
wavelengths where the SF is dominant. The adopted procedure is 
described below. 



3.1 Star Forming Regions in the Host Galaxy 

In order to properly account for the SF emission we need SF 
templates that represent a large range of properties of SF galax- 
ies. We used Spitzer/IRS spectra of 5 different SF galaxies: M82, 
UGC 2982, NGC 6090, IRAS 22491-1808, and Arp 220. These 
galaxies span a wide range of IR luminosity (Lm^^ IO'^'^^'^'^Lq) 
and PAH luminosity (e.g., LpAui.Tjim ~ 

j08.4-9.8^^) jyjQ^j impor- 
tant, we need to take into account the large range and large scatter in 
the ratio between the FIR and MIR lu minosities, known to exist in 
SF galaxies (e.g., ISanders et al1l2003h . As can be seen in Figure[T] 
the scatter in /60//25 (the flux ratio) is between ~3-13. This scatter 
introduces an uncertainty in the determination of the intrinsic AGN 
SED for a specific object (see further discussion in ^4. lb . For each 
of our 5 SF galaxies we choose 8 values within this range (3-13) 
and simulate a 60 flm data point accordingly. We use the simu- 
lated 60 iim data points to construct, from each MIR spectrum, 8 
SF templates covering the NIR to FIR wavelength range. In total. 



Figure 2. Spitzer/IRS spectra of 5 SF galaxies that were used to construct 
the SF templates. Each spectrum is supplemented with 8 simulated data 
points at 60 ^m to sample the observed range in /60//25 shown in Fig.[T] 



our set of SF templates consists of 40 (5 x 8) templates covering a 
large range of IR luminosity, FIR/MIR ratio, and PAH intensities. 
Finally, the observed spectra of the templates are smoothed to avoid 
any systematic features in the subtracted AGN spectra. The 40 SF 
templates are shown in Figure|2] 

We subtract each template from the observed NIR-NIR SEDs 
using the intensity of the PAH emission features. The PAH emis- 
sion is assumed to be solely due to SF in the host galaxy. The 7.7- 
flm PAH feature was detected in 50 sources in our sample. This 
feature is strong and can be used as a reliable indicator for SF activ- 
ity. Recent studies have shown that the ratio between the 11.3 and 
7.7-/im features in A GN may differ significantly from those ob- 
serve d in SF galaxies l lDiamond-Stanic & Riekdr201(ll : ISales et al] 
I2OIOI) . To further test this issue we compared the luminosity ratio 
of the 11.3 and IJ-fim features in our sample to our SF templates. 
All our sources exhibit a much larger ratio than observed in the 
templates, suggesting that either the 1 1.3 flm feature is stronger in 
AGN compared with SF galaxies at similar IR luminosities, or that 
the AGN activity may suppress the shorter wavelength feature. The 
measurement of both features is problematic. The 1.1- flm feature is 
often blended with the adjacent 8.6-/im feature, and the 11.3-/xm 
feature is strongly affected by the adjacent W-flm broad silicate 
feature. 

For each SF template we choose a normalization factor and 
require that after the subtraction the flux at the wavelengths corre- 
sponding to the PAHs will not exceed the noise level. Application 
of this criterion to the 11.3-/im feature resulted in deep "absorp- 
tion" features in the subtracted spectra at the wavelengths of all the 
other PAH features (6.2, 7.7, and 8.6-/xm). Applying the criterion 
to the 7.7-^m feature resulted in residual flux only in the 1 1.3-/im 
wavelength range. We therefore choose to apply this criterion to the 
l.l-flm feature. This can contribute to the uncertainty in the deter- 
mination of the SF contribution and the shape of the intrinsic AGN 
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SED. We further discuss this issue in 94. 1 I For objects with no PAH 
detection the subtraction is limited according to the level of noise 
in the relevant wavelength ranges. 

Given the large scatter in /(,() / fi^ , subtraction of different tem- 
plates using the same criterion at short wavelengths may imply very 
different contribution of the SF component to the FIR emission. 
Given our assumption that the FIR emission in AGN is dominated 
by SF in the host galaxy, we prefer subtractions in which the SF 
contribution to the 60-/lm flux is above 80%. The combination of 
all templates that satisfy both criteria at short and long wavelengths 
define the acceptable range of SEDs. As explained in ^ there is 
another source of uncertainty related to the low spatial resolution of 
the IRAS observations. This uncertainty propagates into the deter- 
mination of the SF component and as a consequence affects the de- 
termination of the intrinsic AGN SED (see 34. H and |4.4| for more 
discussion). 

We calculate the SF contribution to the MIR luminosity by in- 
tegrating over the scaled SF template between 2-35 flm. The me- 
dian contribution is 20% with a wide distribution between and 
75%. 15 sources in our sample exhibit a SF dominated MIR SED, 
i.e. 10 flm silicate feat ure in absorption and very strong PAH emis- 
sion features (see also iDeo et al l |2009l) . Indeed, for these sources 
the SF contribution is much higher compared with the rest of the 
sample. Moreover, SF dominated SEDs are more common among 
NLSIs (12 sources) compared with BLSls (3 sources). This is con- 
sistent with the findings of SIO that NLSIs often exhibit much 
stronger SF activity compared with BLSls for the same AGN lu- 
minosity. SF dominated MIR SED may also be the result of weak 
AGN contribution to this region in the SED due to a small covering 
factor We return to this issue in 34.51 

Subtraction of SF templates with large MIR contribution (> 
35%) is problematic. In these cases the SF contribution, in certain 
narrower wavelengths ranges (e.g. 6.5-7.5 /xm), can exceed the ob- 
served flux, thus the remaining AGN SED becomes meaningless. 
The fact that the SF template can exceed the observed flux at cer- 
tain wavelengths while the integrated flux over the entire 2-35 flm 
range can be as low as 35%, is related to the different SED shapes 
and PAH intensities in SF galaxies. We omit these sources from 
the following analysis of the intrinsic AGN SED. The subtraction 
of the SF template has profound effect at wavelengths longer than 
~ 5-flm, and very little effect on the hot dust component. 

3.2 Clumpy Torus with Silicate-type Dust 

The first AGN-related component is a dusty torus surrounding 
the central en e rgy so urce. We use the clumpy torus models of 
iNenkova et alj j2008 a^ in a similar way to that described in M09. 
The model accepts Lboi as an input to calculate the normalization 
of the SED. Unfortunately, the normalizations that were calculated 
before September 2010, and used in M09, were found t o be off 
by a factor of ~2-3 (see erratum bv lNenkova et aLlboiOl) . By re- 
peating exactly the M09 fitting procedure for the current sample, 
and using the corrected normalization factors, we could not fit the 
spectra with low inclination angle models, i.e. those representing 
type-I AGN. Such corrected models emit more flux than actually 
observed. This suggests that in these cases the torus is exposed to 
too much incoming radiation. 

A possible solution to the above problem is related to the hot 
dust component found in M09. Such a component, located between 
the central source and the inner edge of the torus, must affect the ra- 
diation transfer and the energy balance in the system. This compo- 
nent can have a relatively large covering factor causing a significant 



fraction of the incoming radiation to be absorbed and re-emitted at 
NIR wavelengths. This component is different from the "standard" 
clumpy torus clouds since the temperature of its dust is too high for 
silicate grains to survive (see 94. 2t . Thus, the clumpy torus model 
must be amended to include a pure-graphite dust component in or- 
der to provide a full solution to the observed spectrum. 

To overcome the difficulty associated with the incomplete 
torus model we added an additional free parameter to the fitting 
procedure which in practice is a free normalization parameter This 
parameter represents the fraction of the total radiation reaching the 
clumpy torus i.e., the fraction that is not absorbed by the hot pure- 
graphite dust component. The new procedure does not provide a 
full solution since it does not include the different SED shapes 
"seen" by the torus, but only accounts for the total energy. It pro- 
vides an effective way to determine the contribution of the torus 
component to the IR SED (the torus covering factor). 

The main parame ters and assum ptions of the clumpy torus 
model are detailed in iNenkova et al. I ( l2008iJ) and M09. Unfor- 
tunately, the additional free normalization of the torus makes it 
impossible to fully constrain the geometrical parameters of the 
clumpy torus. However, the contribution of the torus to the MIR 
emission can still be reliably determined. The covering factor of 
the torus, CFjgrus, can be deduced from the ratio between its to- 
tal luminosity and the bolometric luminosity of the central source. 
This definition of the covering factor is equivalent to f{i) in M09 
(eq. 5 there). 

3.3 Dusty NLR 

The second component of the model represents a collection of 
du sty NLR clouds. The mo tivation for this component is explained 
in ISchweitzer et al.l ( 1200 St) and in M09 where it was shown that 
such a component can contribute significantly to the MIR flux of 
luminous AGN. The properties assumed here for these clouds are 
similar to the ones used in M09. We assume constant column den- 
sity clouds with Nil ~ 10^''^ cm^^. We further assume constant 
hydrogen density of 10^ cm^^, solar composition, galactic dust-to- 
gas ratio, and ISM depletion. The important physical parameters for 
this component are the cloud-central source distance (which deter- 
mines the dust temperature), the incident SED and the dust column 
density. For more information on this component see M09. 



3.4 Hot Pure-Graphite Dust 

The third component represents a collection of dusty clouds of gas 
located at the inner edge of the torus. The motivation for this com- 
ponent is explained in M09 where it was shown that such a com- 
ponent is necessary t o expl ain the NIR emission of type-I QSOs. 
iMor & Trakhtenbroll ilOUh have shown that a hot dust compo- 
nent is present in > 80% of type-I AGN and is significantly lu- 
minous. The pure-graphite dust must be external to the broad-line 
region (BLR), where dust cannot survive, and internal to the "stan- 
dard" clumpy torus, where the distances are large enough to allow 
silicate-type grains. 

The physical scale of the hot-dust region is set by two different 
conditions. The outer boundary is the sublimation radius appropri- 
ate for a "typical" dust composed of both (average size) silicate and 
graphite grains. This is given by 

Lbol V^VlSOOKf 
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The inner boundary is the sublimation radius of pure-graphite 
grains with a sublimation temperature of 1800 K, 



:0.5 



^bol 



1046 £ 



1/2 



/1800K\ 



2.8 



pc. 



(2) 



'ergs '/ V ^sub 

This is the innermost radius where an averaged-size graphite dust 
grain can survive. Here we suggest that the pure-graphite dust is 
located in BLR clouds that are between the silicate and graphite 
sublimation radii. We note that these radii depend on the assumed 
grain properties, mostly their size, and are not to be taken at face 
value but rather as typical scales of the system. 

We investigate the conditions in this region by calculating pho- 
toionization models that extend all the way from the innermost 
BLR to the dusty torus. The numerical code used for the calculation 
of the IR continuum emitted by the clouds, as well as their emission 
line spectrum, is the 2009 version of the code ION which is de- 
scribed in detail in various earlier publications (e.g.. iNetzeJIlOOd . 
and references therein). The code includes all the relevant atomic 
and dust-related processes and was compared with CLOUDY (by 
G. Ferland) to ensure similar results under a large range of condi- 
tions including dusty and dust-free gas. 

We adopted the "cloud model" of the BLR which we con- 
sider to be the one most relevant to the situations considered here. 
This is one of three possible generic BLR models addressed in the 
literature. It includes a system of clouds that preserve their mass 
as they move in or out. The properties of the clouds depend on 
the radial coordinate r through an external parameter, s, that de- 
scribes the external confining pressure (presu mably magnetic in 
origin). The model was originally suggested bvlRees et al.l 
and was extended to more realistic cases in iNetzed ( ll99ol) and in 
iKaspi & Netzeil ( Il999h . The main model ingredients are the gas 
density in the clouds {N{r) o= r^'), the cloud size (R{r) oc r"/^), the 
column density (N^oiir) °^ r^"!^), and the geometrical cross section 
of the clouds (A(r) R{r)^). An additional parameter, p, specifies 
the volume density of the clouds, inc{r) oc r^P) and thus the radial- 
dependent covering fraction of the system, C(r). Using the above 
definition we get, dC(r) <>= r^^l^^^dr. 

The actual calculations assume constant density clouds, the 
parameters s and p, and the inner boundary conditions on A' and 
VI col - These can be adjusted to give the best agreement with the ob- 
servations of a specific AGN. Kaspi & Netzer ( 1999) investigated 
a large range of dust-free BLR properties in attempt to fit the time- 
variable emission line spectrum of NGC 5548. They found that the 
values that are consistent with the observations are 1 <s< 1.5 and 
\<p<2. 

Two other types of BLR models were not considered in this 
work. The first is the "local optimally-emitting clo uds" model 
(LOG; e.g., 'Baldw in et alj Il995[ ; iKorista et al.l Il997i) . which as- 
sumes a large range of conditions (density, column density etc.) 
at every location in the BLR. The second is the wind model (e.g., 
lMurravetal]|l995l) which is based on the assumption that BLR 
clouds condense out of a fast outflow, which originates very close 
to the central source, perhaps in the vicinity of the central accre- 
tion disk. The emission line characteristics of the wind model have 
never been compared in detail with specific AGN observations. Our 
preference of the cloud model is based in part on several recent ob- 
servations that directly detect cloud-like objects in the general re- 
gion of the BLR in l ow luminosity AGN (see e.g.. lMaiolino et akl 
l20ICl ; lTurner & Miller.2009. and references therein) 

The parameters of the generic model which has been used in 
this study are: .s = 1, p = l,N(r = Rj c) = 10^'^ cm"-^ andA'„,/(r = 




1 10 
wavelength [urn] 



Figure 3. Normalized spectra, per unit covering factor, of the hot pure- 
graphite dust component. The gas composition is 2 Z; and the cloud prop- 
erties ai'e given in the text. All the spectra shown here are calculated for 
the case of L5100 = 10'**ergs^'. The smallest and largest cloud-source dis- 
tances in the grid are marked next to the lines. The dotted line is the mean 
spectrum of all dusty clouds weighted by their covering factor 



^d.c) = 10^^'66 cm ^, where R^/ c is the graphite sublimation ra- 
dius (eq. |2ll. This combination gives U{r = Rj.c) = 1.6 x 10^^ 
where U is the ionization parameter. The model further assumes gas 
composition of 2 Z,7, and ISM-type depletion. The incoming SED 
is typical of high luminosity AGN and is made of a combination of 
a central accretion disk and a high energy X-ray power-law with a 
break at 50 keV. The results of the calculations depend mostly on 
the value of Li,oi and not on the exact SED shape. 

While the local emission of the graphite dust depends only on 
the grain properties, the dust temperature inside the cloud varies 
by a large factor because of the local grain opacity. Thus a single 
cloud spectrum looks like a combination of many modified black- 
bodies. The calculations provide a grid of dust-produced SEDs that 
are used in the fitting of the hot dust component. Such a grid is 
shown in figure [3] The calculations also provide the emission line 
spectrum of the ensemble of dusty clouds (see 94.2.21 1. 



3.5 Model Fitting and Fit Quality 

The fitting procedure starts with the SF template subtracted SED 
that is assumed to represent the intrinsic AGN continuum. Al- 
though the SF component is subtracted prior to the model fitting, 
it introduces another degree of freedom to the procedure. To fit the 
SEDs with a 3-component models we use minimization in a 
similar way to that described in M09. The only difference is that 
the normalization of the torus component is now a free parameter. 

The fitting algorithm computes values for all possible com- 
binations of torus, NLR and hot dust models. There are seven free 
parameters in the torus model, the normalization parameter and six 
others that describe its geometrical properties (§ [3j- In the NLR 
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Figure 4. Best fit using the three-component model for two representative cases, the NLSl Mrk 896 (left diagram) and the BLSl Mrk 1146 (right diagram). 
Top panels show the best fit model (magenta) and the observed (grey) and subtracted and binned (black) spectra. We also show individual components: torus 
(green), NLR (red) and hot pure-graphite dust (blue). In the bottom panels of each diagram we show the quality of the fit in each wavelength bin, by calculating 
the ratio between the model and the fitted data. 



model there are two free parameters, the cloud distance and a nor- 
malization factor. The distance is changed in steps of 0.075 dex 
between 1 and 850 pc for a source with Lboj = 10*^ erg s^'. For 
the hot dust, there are also two free parameters similar to the NLR 
component, the distance to the cloud and a normalization factor. 
For example, for a source with L5100 = 10^''ergs^' the range of 
distances is changed in steps of 0.03 dex between 0.4 and 6 pc. 

Fig. |4] shows the best fit models for two representative cases, 
the NLSl Mrk 896 and the BLSl Mrk 1 146. The top panel of each 
diagram shows the best fit model (magenta) and the observed (grey) 
and subtracted and binned (black) spectra. We also show individual 
components: torus (green), NLR (red) and hot pure-graphite dust 
(blue). In the bottom panels of each diagram we show the quality 
of the fit in each wavelength bin, by calculating the ratio between 
the model and the fitted data. 



4 RESULTS AND DISCUSSION 

The sample described here was collected from the Spitzer archive 
with no specific selection criteria. Thus, some sources lack certain 
data and/or sufficient quality to be proper analysed. 78 sources in 
the sample (38 NLSls and 40 BLSls) have NIR photometry. In 
the following analysis related to the hot dust component ( g4.2l and 
g4.5l l we only consider these sources. The IRS spectra of 15 sources 
in our sample have very low signal to noise ratio. All silicate and 
R^H features in these sources are undetected and it is impossible 
to constrain the continuum shape. These sources are not included 
the analysis of the intrinsic AGN SED ( PH , the torus SED ( ^431 l, 
and the NLR properties { ^AA\ . 



4.1 Star Formation in the Host Galaxy and the Intrinsic 

AGN SED 

Our method to construct the SF templates is driven by the intrinsic 
scatter in the ratio between the FIR and MIR fluxes found in SF 
galaxies. The subtraction procedure uses two criteria to determine 
the best SF template and its normalization. The first is the condi- 
tion that after the subtraction, any remaining flux at the wavelength 
range of the PAH features would be consistent with the noise level. 
Since our SF library consists of templates that span a large range of 
L/s, only small adjustment of the normalization is needed, and the 
normalization values are usually close to one. The second criterion 
is that the SF template would dominate the FIR part (60-/lm) of the 
SED. 

The SF luminosity (Lsf) of each source is calculated from 
the subtraction process using the integrated IR luminosity of the 
template and its normalization. Figure [5] shows Lgf versus Lboi 
for our sample with different symbols for NLSls and BLSls. We 
also show the Lsf = ^bol line and the relationship Lsf = 10'*"' x 
( Lhni/lO^ ' ^) '^e rgs~' which is adopted from the correlation shown 
in lNetzeJ ( l2009h . 

Altogether, 85 sources (37 NLSls and 48 BLSls) have high 
S/N IRS spectra and AGN dominated MIR SEDs. For each of these, 
the range of different templates that meet the PAH criterion is rel- 
atively narrow, and the uncertainty on the shape of the AGN SED 
is small. The SF-template subtracted spectra are normalized at 14- 
jlm. and then used for constructing median intrinsic AGN SEDs 
for two groups divided by FWHM(Hj8), i.e. NLSls and BLSls. 
The choice of the normalization point may have some effect on the 
shape of the median SED. We choose to normalize the spectra at 
14-/im since this wavelength is less effected by the broad silicate 
emission features and the large scatter at longer wavelengths. The 
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Figure 5. Star fomiation luminosity in AGN liost galaxies vs. AGN 
bolometric luminosity for the two groups of NLSls (circles) and BLSls 
(squares). The SF luminosity is calculated using the IR luminosity of 
the chosen template and its normalization. The solid line represents 
LsF = Lbol and the dashed line represent the relation Lgp = 10*^ x 
(Lbol/10'")'"£rgs ' adopted from the correlation shown in lNetzeJ J2009l) . 



Table 1. Median Intrinsic AGN SED 



Rest Wavelength 
(Mm) 


Intrinsic AGN SED 
AL;[ (arb.units) 


0.51 


0.8401 


1.2 


0.9956 


1.7 


1.0182 


2.2 


1.0364 


2.7 


1.0479 


3.2 


1.0555 



Table[T]is published in its entirety in the electronic edition 



right panel of Figure |6] shows the median AGN SED for the two 
populations of narrow- and broad-hne AGN in our sample. The left 
panel shows the median AGN SED for all 85 sources. Examining 
Figure [6] it is evident that the general shape of the SEDs, regard- 
less of line widths, is very similar. All SEDs exhibit three distinct 
peaks at short (< 5-/im) and medium (-^ 10 and 20-/xm) wave- 
lengths. The 10 and 20-/im peaks are likely dominated by silicate- 
dust emission. The median intrinsic AGN SED for all 85 sources 
is also listed in Table [T] We also examined the median AGN SEDs 
for two luminosity groups defined by the median luminosity of the 
sample at 6-/im (lO'*'*'^ erg s^'). We do not find any significant 

difference between the two SEDs. 

Our result is consistent with the shape found bv lNetzer et al.l 

( l2007h for the QUEST sample. The QUEST sample represents the 
high end of the AGN luminosity rang e of our current sam ple, and 
mostly consists of broad line sources. iNetzer et alj ( |2007|) also di- 
vided the QUEST sample into two groups of strong and weak FIR 



emitter s. Again, the shap es of the AGN SEDs were remarkably 
similar. iDeoetalJ ^20091) employed a similar subtraction method 
to a large sample of type 1 and 2 Seyfert galaxies. These au- 
thors find that the AGN continuum emission drops rapidly beyond 
^ 20-Atm for a l l AGN types, regardless of SF activity. Recently, 
iMullanev et al.l ( 1201 1[) employed a somewhat different method to 
determine the intrinsic AGN SED of a medium-size sample of type- 
I and type-II AGN. These authors fit simultaneously the IR SEDs 
using a SF template and a broken power law that represents the 
AGN emission. They find that the warm dust peak is located at 
somewhat longer wavelengths, between ~ 15 and 45-/im, before 
dropping rapidly towards FIR wavelengths. This result is not con- 
firmed by our analysis probably due to the different way of ac- 
counting for the SF contribution, and the different nature of the 
MIR ingredients they consider (a power-law). We suspect that it 
may also be affected by the presence of type-II AGN in their sam- 
ple and the difficulty to properly account for the host contribution 
in these sources. 



4.2 Hot Dust Properties 

A major goal of this work is to explain the NIR emission in type-I 
AGN and to identify the physical properties of the component re- 
sponsible for this emission. The hot dust models used here provide 
both the continuum emission, which peaks at NIR wavelengths, and 
the UV-NIR lines emitted by this component. 

4.2.1 Continuum emission and distances 

The most important parameter of the hot dust component deter- 
mined by the fitting procedure is its luminosity (Lhd)- Since the hot 
dust is optically thick at all UV-optical wavelengths, Lhd is simply 
a measure of the covering factor of this component, CFhd- Figure|7] 
shows a strong correlation between Lhd and L(,oi- NLSls (blue cir- 
cles) exhibit a similar trend to that of BLS Is (red circles) although 
the scatter in this group is larger. The solid line in Fig. |7] repre- 
sents a simple least squares fit to the data of t he current work with 
a slope of 0.9. Fig.|7]also shows the results of iMor & TrakhtenbrotI 
f201 1) (grey dots). These authors fitted the (observed) optical to 
MIR SEDs of ~ 15000 high luminosity QSOs and measured the lu- 
minosity of the hot dust component. The relation between Lhd and 
Lboi, over a large luminosity range, indicates that CFhd in type-I 
AGN spans a relatively limited range (~0. 1-0.2). The NLSls that 
lie much below the relation of Fig. |7] may represent sources with 
exceptionally small CFhd (see g4.5t 

The second parameter that is determined by the fitting pro- 
cedure is the distance from the centre to the hot dust compo- 
nent (i^HD)- Figure [8] shows this distance against L\,q\. The sil- 
icate and pure-graphite dust sublimation radii defined by equa- 
tions [T] and (2] respectively, are also shown. We also show the re- 
lation between the emissivity-weighted BLR radius for the Hj3 
line, i?BLR(Hj8) an d Lhni- This radiu s is based on an up-to- 
date version of the Kaspi et al. ( 2005) relation taking into ac- 
count the modifications of Bentz et al. (2009): ^blr(H/3) = 0.35 x 
(L5100/IO'** ergs^')"'^^ pc. Fig.[8]demonstrates that the clouds are 
clearly situated outside the dust-free BLR and inside the edge of the 
"standard" silicate-dust torus. 

Several other model parameters affect the derived distance, for 
example the column density of the graphite dust (which is related to 
the gas metallicity) that affect the mean dust temperature in an op- 
tically thick clouds. However, changing the values of these param- 
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Figure 6. Intrinsic AGN SEDs. Right: median intrinsic AGN SEDs for two sub-groups of NLSls (blue) and BLSls (red). Dashed lines represent the 25th and 
75th percentiles of each sub-group. Left: median intiinsic AGN SED for all 85 sources (see data in Table[T). As in the right panel, dashed hnes represent the 
25th and 75th percentiles. All SEDs have very similar shape and exhibit three peaks at short (< 5-^m) and medium (~ 10 and ~ 20-^m) wavelengths. The 
SEDs are normalized at 14 /im and their shape may be somewhat effected by the chosen normalization point. 
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Figure 7. Luminosity of the hot dust component (Lhd) vs. AGN bolometric 
luminosity (Lbol)- Results from the current work are shown as circles and 
squares (NLSls and BLSls, res pectively). Results for a much larger sample 
of high luminosity QSOs from iMor & Trakhtenbroll j201lh are shown in 
grey. The black line represents a simple least squares fit to the data of the 
current work. The 37 sources without NIR data are not included in this plot 



eters do not significantly alter the results found here. For example, 
assuming solar metallicity, instead of 2 Z0, enlarges the derived 
distances, for all sources, by about a factor of 1.25. These distances 
are still consistent with the assumed location of the dusty clouds. 



4.2.2 The emission line spectrum of dusty BLR clouds 

So far we have only discussed the hot dust continuum emission. 
However, these "dusty-graphite-BLR" clouds also produce broad 
emission lines that will be observed in the UV-NIR part of the 
spectrum. This aspect has never been investigated and is likely 
to result in more observational constraints on the hot dust model. 
Fig. [9] shows normalized line fluxes, per unit covering factor, for 
the 2 Zq) metallicity BLR model considered here. The line inten- 
sities are shown as a function of distance for the specific case of 
L5100 = 10'*^ergs^' . Other cases can be obtained from the diagram 
by scaling r with L'/^ 

Given the known radial changes in the differential covering 
fraction, dC(r), we can integrate over the line emissivity to ob- 
tain the cumulative line luminosities produced by the dusty-BLR 
clouds. This is shown in Fig.[TO] The diagram illustrates the large 
increase in line flux up to the graphite sublimation radius and 
the much slower increase beyond this boundary. For the model in 
question, the covering fraction of the dusty-graphite-BLR clouds 
is almost exactly the same as the cumulative covering fraction of 
the dust-free clouds. If the dusty-graphite-BLR clouds were dust 
free, the total emitted flux of lines such as Lya and Hj8, would 
roughly double over this region. Since much of the ionizing radi- 
ation is absorbed by the graphite grains, the increase in intensity 
is significantly reduced. For example, in the dusty-graphite-BLR 
model considered here, the increase in Lya is 11%, in H/3 7%, in 
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Figure 8. Distance to the hot pure-graphite dust components (Shd) vs. 
AGN bolometric luminosity (L|,ol) for two sub-groups of NLSls (circles) 
and BLSls (squares). Dashed and solid black lines represent the sublima- 
tion radii of silicate and pure-graphite dust (eqn.[T]and[2), respectively. Dot- 
dashed line represent the fiBLR(H/3)-Lboi relation (see text). 
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Figure 9. Normalized line emissivity per unit covering fraction for strong 
broad emission lines for the case of twice solar metallicity. The graphite 
sublimation radius in this case is at 1.32 pc and the silicate sublimation 
radius at 3.47 pc. These distances are clearly visible due to the large drop 
in the emissivity of all lines resulting from the absoiption of the ionizing 
radiation by the grains and the selective depletion of most metals. 
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Figure 10. Cumulative emission line intensities as a function of distance 
from the central source for the model presented in Fig. |9] The dashed ver- 
tical lines represent the graphite and silicate dust subhmation radii for the 
model. For most lines, the increase in total emitted flux is very small al- 
though not completely zero (see text for details). 

C IV 11549 6%, in Mg II 12798 38%, iii He 1 1 10830 69%, and in 
Fe II only 3%. Thus, the only significant contributions to the total 
line intensities are due to Mg II and He I lines. 

The BLR model presented here is not entirely consistent with 
the SED fitting procedure of ^since all clouds in the pure-graphite 
zone are contributing to the emission lines yet the SED fitting takes 
into account only one such cloud. As illustrated in Fig. [3] this is not 
a severe limitation since the covering-factor-weighted mean spec- 
trum has an SED that is very similar to that produced by a single 
cloud. 

The BLR models considered here do not represent the entire 
range of possible properties. In particular, the effect of the internal 
dust on the line emission depends on the ionization parameter ((7) 
since at lower V , the fraction of the ionizing continuum absorbed 
by the dust is reduced. This can result in a larger increase of the line 
intensity in this zone. The exact values of f/, the gas density, and 
column density of the cloud model are determined by the parame- 
ter s ( g3.4t . A more detailed consideration of these possibilities is 
deferred to a future publication. 



4.3 Torus Properties 

As explained in 33.21 the presence of hot dust clouds at the inner 
edge of the torus will affect the torus emission spectrum. The hot 
dust clouds will absorb a large fraction of the incident UV radiation 
and re-emit it at NIR wavelengths and at different directions than 
the original AGN radiation. This will change the SED and total 
luminosity incident on the torus. 

As a first attempt to solve this issue without resorting to a full 
new torus model (which is not yet available), we subtracted the 
hot-dust luminosity, Lhd> from Lboi and used this value in order 
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to calculate the torus models (see 93.21 ). For this reduced amount 
of incident luminosity, we used the exact same scheme as M09 in 
which the torus normalization is not a free parameter but rather set 
by the (now reduced) AGN luminosity. We found that most sources 
still could not be fitted using low inclination angle, those represent- 
ing type-I AGN tori. To fit the spectra using such low inclination 
models one would need to artificially reduce Lbol by a factor that 
is larger than the one implied by the hot dust emission. This unsuc- 
cessful attempt may imply that the viewing angle of the hot dust 
(which is not considered here) can play an important role or that 
the geometry of the hot dust clouds is different from that of the 
torus. For example, if the hot dust component has similar geometry 
to that of the torus, its emission would also be highly anisotropic. 
Moreover, if the two components have different inclination angles, 
the amount of radiation that the torus would receive from the hot 
dust component would be different than the amount inferred by the 
observed NIR flux. We conclude that the only way to fully examine 
this possibility is to include a pure-graphite dust component into 
the radiative transfer calculation of the clumpy torus. 

In further attempt to solve this issue, we set the normalization 
of the torus component to be a free parameter of the fitting proce- 
dure ( 93. 2t . This way, models with low energy output, small phys- 
ical size, or high inclination angles w.r.t. the line of sight, fit well 
the observed spectra. This adds a considerable uncertainty to the 
derived values of the geometrical parameters of the fit. However, 
the average sample properties are less likely to be affected. 

The following conclusions refer to the distribution of the dif- 
ferent torus parameters for all sources with good IRS S/N ratio (49 
NLSIs and 51 BLSls). We found that some parameters have a nar- 
row distribution around a mean value and in some cases only a sin- 
gle acceptable value. Other parameters exhibit a broader, more uni- 
form distribution. The parameter distributions of all sources were 
combined by giving each value within the acceptable range in a 
certain source its relative weight in the distribution. The results are 
(see M09 for more explanation of the parameters): 



(i) The median torus width parameter is a = 35. 

(ii) The radial extent of the torus, Y, has a broad distribution 
with a median value of 30. The range in Y implies torus outer sizes 
that range between --^1.7 to 85 pc. 

(iii) The average number of clouds along an equatorial ray, A'o, 
has a broad distribution with a median of 4 clouds. 

(iv) The parameter that specifies the radial power-law distribu- 
tion of the clouds, q, has a median value of 1.7. This parameter is 
related to the anisotropy of the torus radiation. As q decreases the 
torus radiation becomes less isotropic. 

(v) The median optical depth of a cloud is Ti/=38. This parame- 
ter has a broad distribution with Xy > 20. Since the requirement for 
large MIR optical depth is Ti/~ 10, it is not surprising that the torus 
models are not very sensitive to the exact value of Ty beyond this 
value. 

(vi) The torus inclination angle shows a broad distribution for 
all ; < 60° with a median of 50 ° . This again is consistent with the 
assumption that the direct line-of-sight to the centre of type-I AGN 
is almost completely free of obscuring material. 



We stress again that we do not attach great significance to a specific 
value of a specific torus parameter but rather make sure that all 
parameters lie within a reasonable range for type-I AGN. 
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Figure 11. NLR cloud distances vs. AGN bolometric luminosity. NLSIs 
are marked by circles while BLSl are marked by squares. The slope from 
M09 (0.47) is shown by the grey line for comparison. The 15 sources with 
very low IRS S/N ai'e not included in this plot. 



4.4 NLR Properties 

The important parameters of the NLR component are the distance 
from the centre, which determines the dust temperature, and the to- 
tal dust luminosity which determines the fraction of MIR flux emit- 
ted by such clouds. The properties of this component are similar 
to those described in M09. Figure [TT] shows NLR cloud distances 
against the bolometric luminosity of the objects. Using simple least 
squares regression we find the following scaling relationship. 



:460x 



^bol 



0.67±0.05 



pc. 



(3) 



10^6 erg s- 

The grey line in Fig.[TT]represents the scaling relation (slope 0.47) 
found for the QUEST sample by M09. Although we find a slightly 
steeper slope, both results are consistent within the observed scat- 
ter. The distances found here are similar to the values found by 
M09, however the scatter in the current sample is larger. 

As explained, the main difference from M09 regarding the 
NLR component is the more detailed treatment of the SF contri- 
bution. For example, the emission from the dusty NLR component 
peaks at MIR wavelengths (above ~20-/im) where the SF contribu- 
tion to the observed spectrum is substantial. M09 used a single M82 
template to subtract the observed spectrum. For sources with strong 
SF activity (i.e., those with high Lboi or NLSIs that exhibit strong 
SF, see SIO) this template does not properly account for the ob- 
served FIR emission. Other SF templates have larger contribution 
towards shorter wavelengths (~30-/im). Consequently, the NLR 
component would have less weight in the AGN MIR spectrum and 
different spectral shape. This can result in larger NLR distances and 
smaller luminosity. We consider this as an additional uncertainty on 
the determination of the NLR distance and its relative contribution. 



Graphite Dust in AGN 1 1 




0.2 0.4 0.6 0.8 1 

Covering Factor 



Figure 12. Cumulative distribution functions of tlie covering factors of tlie 
three AGN dusty components. For all components, the covering factors tend 
to be smaller in NLSls (blue) than in BLSl(red). The distribution of CFhd 
and the torus covering factor are very similar indicating that about half of 
the obscuration is due to hot pure-graphite dust clouds which are typically 
not included in "standard" torus models. Dashed lines in the upper panel 
represent CFhd obtained by fitting the 1-5 ;tlm range using the hot dust 
component only (without a torus component - see text) 



4.5 Covering Factors of the Dusty Components 

A major assumption of this work is that the entire MIR spectrum, 
after starburst subtraction, is reprocessed AGN radiation. This can 
be used to deduce the covering factor of the central source by the 
three components (see also M09 and references therein). The cov- 
ering factors are defined by the ratio between the total luminosity of 
each component and Lboi- Figure[T2]shows the cumulative distribu- 
tion functions of the covering factors of the different components. 
The blue line represents the NLSls in the sample and the red line 
represents the BLSls. As seen from the diagram, NLSls tend to 
have smaller covering factors compared with BLSls. The median 
value of CFhd for the NLSls is 0.23, for the BLSls 0.27, and for 
the entire sample is 0.25. The median value of the torus covering 
factor is 0.24 for the NLSls and 0.33 for the BLSls. The covering 
factor of the NLR component has a median value of 0.03 for the 
NLSls and 0.07 for the BLSls. 

SIO found that NLSls tend to have stronger SF activity com- 
pared with BLSls at a given AGN luminosity. Since we fit the SF 
subtracted SED this may effect the derived torus and NLR covering 
factors. However, this cannot explain the difference in the covering 
factors of the hot dust component since the SF contribution to the 
1-5 jlm wavelength range is very small. Host contribution to the 
optical range (i.e. to L5100) would result in an overestimation of 
the AGN bolometric luminosity and consequently lower the esti- 
mated covering factor (see discussion on the uncertainty of Lboj 
in Such host contribution to the optical spectrum however, is 
more significant in lower luminosity AGN regardless of their line 
widths, and cannot explain the observed discrepancy. 



We further checked the influence of the torus component on 
the fitting of the NIR wavelength range (1-5 /im) where the hot 
dust component dominates the spectrum. We fitted this limited 
wavelength range using a hot-dust component only, without in- 
cluding a torus component. The dashed lines in the upper panel 
of Fig.[T2]represent the CFhd obtained in this manner. The median 
value of this CFhd for all sources is 0.29 and is slightly higher than 
that found by the fitting of the entire SED using all three compo- 
nents (0.25). Thus using NIR photometry only, CFhd (and Lhd) 
are overestimated, on average, by 16%. 

In order to explore the relation with AGN properties we fo- 
cus on CFhd and compare it against Lboh ^bh. and L/LEdd- 
Fig. 1 131 shows CFhd against these properties for the current sam- 
ple of NLSls (blue c i rcles) and BLSls (red circles). The results of 
iMor & TrakhtenbroJ jioTU grey dots) are shown for comparison. 
The black dashed lines in Fig. [T3] represent the 99th percentiles 
bound aries of the CFhd distribution of the iMor & T rakhtenbro^ 
( I2OI ih sample. These percentiles were calculated by assuming that 
the CFhd distribution, in a certain luminosity bin, should be sym- 
metrical around the peak value and by mirroring the high-CFnD 
side of the distribution. This is done in each 0.2 dex luminosity bin 
separately. Mor & Trakhtenbrot (2011) suggest that all the points 
which lie below the lower dashed lines can be regarded as hot-dust- 
poor AGN. The fraction of such sources in their sample is ^20% 
and does not depend on source luminosity. 

The left panel of Fig. [T3] presents a clear anti-correlation, in 
the sense of a decreasing CFhd with increasing Lboi- This trend is 
further confirmed by both Pearson's and Spearman's rank correla- 
tion tests (p value <C 0.01 for both tests). Severa l earlier studies 
sugge sted a similar trend (e.g., IWang etal.ll2003 :l Maiolino et al.l 
|2007; Treister et al. 2008), however these were based on the total 
MIR emission which translates to the cove ring factor of th e entire 
dusty structure, and not just the hot dust. ICallagher et al.l l l2007h 
suggest that this CFno-ibol anti-correlation may be the manifes- 
tation of dust extinction in the host galaxy Pl lMor & TrakhtenbrotI 
( I2OI l[) . however, showed that such an anti-correlation is still very 
significant even for the bluest sources in their sample. These "blue" 
sources have optical slopes greater than —0.4, and are less likely to 
be affected by extinction. The physical mechanism responsible for 
the decrease of covering factor with L[,oi is st ill undetermin ed. One 
possibility is a "receding torus" scenario ( Lawr enc^l 1 99 ih . where 
higher luminosity implies larger dust sublimation distance, and 
hence an obscuring structure that is located farther away from the 
centre. In this scenario, the geometry of the hot-dust clouds must 
be toroidal and have a constant height. Another possibility is that 
the CEnD-i-bol anti-correlation is analogous to the anti-correlation 
found between the equivalent width of different BLR lines (Hj3 
and C IV) and AGN lumino sity (e.g., lBaldwin|[T977l; I Netzer et al.l 
I2OO4I : iBaskin & Laod[2005h . The CFnD-ibol anti-correlation can- 
not explain the lower CFhd of the NLSls in the current sample, 
since both BLSls and NLSls sub-groups span a similar luminosity 
range. 

The NLS 1 s in our sample have lower Mgn than the BLS 1 s. As 
can be seen in the middle panel of Fig. [13] most of the sources with 
log(MBH/M3) < 7.3 are NLSls and ~20% of them have CFhd 
lower than 0. 1. No significant correlation was found between CFhd 
and Mbh- If the NLSls that have very low CFhd are truly hot dust 
poor AGN it may suggest that these sources are more common in 



' iGallagher et al.l feOOTh present a somewhat different covering factor since 
they integrate over a larger range of IR wavelengths (~1-100 jlm). 
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Figure 13. Covering factor of the hot pure-graphite dust c omponent against l^gi (left ), Mbh (middle), and L/L^^ (left panel). In all panels, NLSls are 
represented by the circles and BLSls by squares. Results of Mor & Trakhte nbroll flOllh (grey dots) are shown for comparison. Dashed black lines represent 
the peak of the CFhd distribution and the 99th percentiles boundaries - see text 



low Mbh sources. iMor & Trakhtenbroll l lioTlh did not find such a 
dependence of the fraction of hot dust poor AGN on Mbh- How- 
ever, the reason that this dependence was not found can be related 
to the relatively narrow Mbh range that the iMor & Trakhtenbroll 
l|) sample span. 

L/Ledd is known to be higher in NLSls. The right panel of 
Fig. [Ts] suggests an anti-correlation between CFhd and L/Lgj^, 
over two orders of magnitude in L/L^^i. The NLSls that lie be- 
low the CFno-ibol and CFhq-Mbh relations a re als o the sources 
with the highest L /Le^h. iMor & Trakhtenbroll ( 1201 ll) did not find 
significant correlation between CFhd and L/LEdd- However, this is 
probably due to the relatively narrow range of L/Ledd in their sam- 

S ile. S imila r to the relation with t he AGN luminosity, iNetzer et al.l 
20041) and lBaskin&Laoil ^20051) found that the equivalent width 
of the broad Hj3 and C IV lines also anti-correlates with L/L^i^. 



5 CONCLUSIONS 

We conducted a detailed investigation of the NIR-FIR SED of a 
large sample of Spitzer-observed AGN. We fitted the spectra of 5 1 
NLSls and 64 BLSls using a three component model: a clumpy 
torus, a dusty NLR, and a hot pure graphite-dust component which 
is a continuation of the BLR. The fitting was performed on SF sub- 
tracted SEDs using SF templates that take into account the entire 
range of possible host galaxy properties. 

The main results of the investigation are: 

• The intrinsic MIR AGN SED is very similar in NLSls and 
BLSls regardless of AGN luminosity. All SEDs exhibit three dis- 
tinct peaks at short (< 5-/xm) and medium (--^ 10 and ~ 20-/xm) 
wavelengths. The first peak is likely dominated by emission from 
pure-graphite dust located outside the dust-free BLR and within 
the "standard" silicate-dust torus. The other two peaks are likely 
dominated by silicate-dust, corresponding to the known 10 and 18- 
/im broad emission features. The long wavelength downturn of the 
AGN SED is at about 20-25 /im. 

• Our detailed modelling of the hot-dust component allows us to 



explore the emission line spectrum of such clouds. Most line emis- 
sion is dramatically suppressed in these dusty clouds. However, the 
Mg II A2798 and He I lines are still contributing significantly to the 
total BLR spectrum. 

• We calculated the covering factors of all AGN-related compo- 
nent and found that the covering factor of the hot-dust component 
is about 0.25. The covering factor of the torus component is com- 
parable with a median value of 0.28. The NLR has a much smaller 
covering factor of about 0.07. The NLSls in our sample tend to 
have smaller covering factors than BLSls. 

• The covering factor of the hot-dust component, CEhd. anti- 
correlates with both Lboi and L/Le(J(J. The CEnD-^bol anti- 
correlation may be related to the idea of a "receding torus" where 
higher luminosity implies larger dust sublimation distance, and 
hence an obscuring structure that is located farther away from the 
centre. This however requires that the hot dust component would 
have a toroidal structure with constant height. The origin of the 
CFno-i/^Edd relation is unknown and is not a consequence of the 
CFHD-i'bol relation since the NLSl and BLSl sub-samples have 
similar distribution of Lboi- Both the CFno-^'bol and the CFhd" 
L/Le(J(J relations may be analogues of the anti-correlations found 
between the equivalent widths of broad H/3 and C IV lines and AGN 
luminosity and L/L^dd- 
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